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Abstract. Cosmological weak lensing gives rise to correlations in the ellipticities of 
faint galaxies. This cosmic shear signal depends upon the matter power spectrum, 
thus providing a means to constrain cosmological parameters. It has recently been 

6 

to. 

with decreasing survey depth. Here we consider the two-point shear correlation 
function, and demonstrate that photometric redshift information can be used to 
suppress the intrinsic signal; at the same time Poisson noise is increased, due 
to a decrease in the effective number of galaxy pairs. The choice to apply such 



proposed that intrinsic alignments arising at the epoch of galaxy formation can 
also contribute significantly to the observed correlations, the amplitude increasing 



a redshift-depending weighting will depend on the characteristics of the survey 
in question. In surveys with (z) ~ 1, although the lensing signal dominates, the 
measurement error bars may soon become smaller than the intrinsic alignment 
signal; hence, in order not to be dominated by systematics, redshift information 
in cosmic shear statistics will become a necessity. We discuss various aspects of 
this. 
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1. Introduction 

The distortion of distant galaxies by the tidal gravitational field of intervening matter 
inhomogeneities has become known as "cosmic shear" . Since this lensing signal depends 
upon the matter power spectrum, it is an important cosmological tool, as was proposed 
in the early 1990s by Blandford et al. (1991), Miralda-Escude (1991) and Kaiser (1992). 
This set the scene for further analytic and numerical work (eg. Kaiser 1998; Schneider et 
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al. 1998; White & Hu 2000), taking into account the non-linear evolution of the power 
spectrum which results in increased power on small scales (Hamilton et al. 1991; Peacock 
& Dodds 1996). 

During 2000, four teams announced the first observational detections of cosmic shear 
(Bacon et al. 2000; Kaiser et al. 2000; van Waerbeke et al. 2000; Wittman et al. 2000; 
Maoli et al. 2001), demonstrating the feasibility of its study. Various statistics are used to 
quantify cosmic shear, and to compare the observations with predictions for cosmological 
models. Here we focus on the shear correlation functions; other measures include the 
aperture mass statistic (Schneider et al. 1998) and shear variance (e.g. Kaiser 1992). 

Besides its dependence on the large-scale structure in the Universe, the magnitude of 
the effect is also sensitive to the redshift distribution of the galaxies used in the analy- 
sis. So far, direct redshift estimates have not been obtained for the samples concerned. 
Typically the mean redshift has been estimated from surveys of similar depth, and a cor- 
responding redshift probability distribution has been assumed (e.g. van Waerbeke et al. 
2001; Hoekstra et al. 2002a). Motivated by the recent interest in obtaining photometric 
redshifts for cosmic shear surveys, we consider how to make use of redshift information, 
and its impact on the constraints permitted by the two-point statistic under considera- 
tion. 

In weak lensing analyses, it is assumed that the background galaxies are randomly 
oriented so that their mean intrinsic ellipticity (e s ) = and any correlation in observed 
cllipticities arises from gravitational lensing. The magnitude of any intrinsic alignment 
of the background source population has recently been the subject of many numeri- 
cal, analytic and observational studies (e.g. Croft & Metzler 2000; Heavens et al. 2000; 
Crittenden et al. 2001 (CrOl); Catelan et al. 2001; Mackey et al. 2002; Brown et al. 2002, 
Hui & Zhang 2002). Due to differences in the assumed origin of intrinsic alignments and 
the type of galaxies considered, the numerical and analytic estimates span a couple of 
orders of magnitude in amplitude. For example, the analytic model of CrOl assumes that 
any intrinsic signal arises from correlations between the angular momenta of galaxies, 
whereas that of Catelan et al. (2001) relies on tidal shear correlations. Nonetheless, all 
studies conclude that intrinsic alignments are more important for shallower surveys (be- 
coming comparable to the lensing signal for a mean survey redshift (z) ~ 0.5) and that 
the correlation falls off quite rapidly with source separation. 

The structure of this paper is as follows. In the next section we outline the rela- 
tionships between the matter power spectrum and the lensing correlation functions. The 
dependence on the redshift distribution and on cosmology is highlighted. Correlation 
functions are derived from observational data, and in Sect. 3 we describe how practical 
estimators are used for this purpose. In Sect. 4 a toy model to account for intrinsic align- 
ments is developed; the amplitude of the correlation is normalised to that of CrOl. The 
introduction of a weighting factor to minimise any contribution of intrinsic alignments to 
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the shear correlation function is considered in Sect. 5. In Sect. 6 we present the results of 
applying such a weighting factor to surveys of different mean redshift. We finish in Sect. 7 
with a discussion of the results and a future perspective of using redshift information in 
cosmic shear analysis. 

For a review of cosmological weak lensing and the relevant aspects of cosmology, 
see Mellier (1999) and Bartelmann & Schneider (2001; hereafter BS01). In addition, the 
present status and outlook for cosmic shear studies is summarised by van Waerbeke et 
al. (2002). 

2. Power spectra and lensing correlation functions 

In this section, the relationships between the power spectra and observable lensing cor- 
relation functions are outlined; throughout, we adopt the notation of BS01. 

The three-dimensional density fluctuation field S(r) at comoving position r is a ho- 
mogeneous and isotropic random field, with a two-point correlation function given by 

(8(r)5(r'))=C 55 (\r-r>\) . (1) 

Such a field is described by its power spectrum Pa(|k|) (k being the comoving wave- 
vector), which is the Fourier transform of the two-point correlation function. In Fourier 
space one can write 

*(k)J*(k')) = (2n) 3 S B (k - k')P tf (|k|) , (2) 



where x* denotes the complex conjugate of x, x denotes its Fourier transform, and Sy> is 
the Dirac delta function. 

The effective convergence R e s(0) depends upon the weighted integral of the density 
contrast along the line of sight 

3# 2 ft m p ^ llrr ^ 6{f(w)0,w] 



= ^r i y o d^H/w^, (3) 

where Hq and f2 m are the values of the Hubble parameter and the density parameter at 
the present epoch, and a(w) is the scale factor at comoving distance to, normalized such 
that a(0) = 1 today. The horizon distance is denoted by wh, and f(w) is the comoving 
angular diameter distance, so that the comoving separation is f(w)9. The function f(w) 
depends on the spatial curvature, K: 

K- 1 / 2 siniK^w) {K>0) 

w (K = 0) ] (4) 

(-X)- 1 /2 sin h[(-i<:) 1 /2 w ] (K<0) 

later we assume that K = (i.e. £l m + £!a = 1). The function W(w) accounts for the 
sources being distributed in redshift, and consequent differences in lensing signal, 

W(w)= I™" dw> p(w') f{w '~ w) =(R(w,w')), (5) 

Jw /K) 



fW = { 
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where p(w')dw' is the comoving distance probability distribution for the sources. The 
function R(w,w') = f(w' — w)/f(w') is the ratio of the angular diameter distance of a 
source at comoving distance w' seen from a distance w, to that seen from w = 0. 

The power spectrum P K (£) of the effective convergence, or equivalently of the shear 
P 1 {£) (see for example BS01), is related to that of the density fluctuations through a 
variant of Limber's equation in Fourier space (Kaiser 1998) 

where i is the angular wave- vector, Fourier space conjugate to 6. 

We see now that P K and P 7 are intimately related to fi m , and to P$ which in turn 
depends on Q\, the shape parameter T and on as, the density fluctuations in spheres 
of radius 8/i _1 Mpc. In addition, note the dependence upon the redshift distribution of 
the sources. Access to the cosmological parameters is provided through the (observable) 
shear correlation functions, which we now turn to. 

The shear has two components and can conveniently be represented as the complex 
quantity 7 = 71 + i72- Throughout, we assume the weak lensing limit so that |— y | <C 1. 
Consider a pair of galaxies at position # and # + 6, where the polar angle of their 
separation vector 9 is </>. The tangential and cross-components of the shear for either of 
these galaxies are 

7t = -He [ 7 e" 2 ^] ; 7 x = -Im [ 7 e" 2 ^] . (7) 

Using the notation of Schneider et al. (2002a), the shear correlation functions are defined 
as 

£±0?) = <7t7t> ± <7x7x> (9) = ^ J 0l4 (#) P K (i) , (8) 

where J„ are n-th order Bessel functions of the first kind. P K (£) can be written in terms 
of the observable correlation functions: 

/•OO 

p k (£) = 2tt dee K ± (fl)Jo, 4 (^)] , (9) 

Jo 

making use of the orthogonality of Bessel functions. From now on we focus on £ + , since 
this contains most cosmological information on the scales of interest. 

Another quantity that is often determined in cosmological weak lensing studies is the 
shear variance inside a circle of radius 9 C , which is also related to the convergence power 
spectrum 

(m>)m = lj«uP M ^. do, 

This can also be determined directly by placing circular apertures onto the data field but 
is affected by any gaps; these may arise, for example, due to the need to mask regions 
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containing bright stars and their diffraction spikes. However, the shear variance can also 
be obtained from £ + : 

(l<H*c) = /^£ + W + (£) , (ID 

where van Waerbeke (2000) showed that 

f i [4 cos" 1 (f ) - aV4-a; 2 | (a; < 2) 
S+{x) = { 71 1 {2) J V " ' . (12) 

[ (x > 2) 

2.1. Choice of cosmology, power spectrum and source redshift distribution 

Unless otherwise stated, our fiducial cosmology is a ACDM model with Q m — 0.3, fl\ = 
0.7 and H = TOkms" 1 Mpc~\ 

A scale- invariant (n = 1, Harrison-ZePdovich) spectrum of primordial fluctuations 
is assumed. Predicting the shear correlation functions requires a model for the redshift 
evolution of the 3-D power spectrum. The fitting formula of Bardeen et al. (1986; BBKS) 
is used for the transfer function, and the Peacock and Dodds (1996) prescription for the 
evolution in the nonlinear regime. The power spectrum normalisation is parameterised 
with & § = 0.9, and the shape parameter T = 0.21. 

As an aside, the differences in predicted shear correlation functions using the fitting 
formula of Eisenstein & Hu (1999) for the transfer function (rather than that of BBKS) 
were found to be minimal. 

The normalised source redshift distribution is parameterised using the form suggested 
by Smail et al. (1995): 

p{z)dz = ^WW) (£) 

where T 7 (a;) is the gamma function. We set (3 = 1.5 so that (z) w 1.5 zq. The value of zq 
is adjusted to give different redshift distributions, which would in practice be obtained 
from surveys with different limiting magnitudes. For example, Fig.|l| shows p(z) for (z) = 
0.5 and (z) — 1.0, roughly corresponding to limiting magnitudes of / ~ 22 and ~ 25 
respectively. 



dz , 



(13) 



2.2. Some dependencies of£ + 

In this subsection we note some dependencies of on cosmological parameters, and on 
the redshift distribution of the sources from which it is measured. 

To illustrate the dependency of £ + on source redshift distribution, Fig.^ shows £ + for 
three values of zq, corresponding to (z) = 0.8, 1.0 and 1.2. As expected, the amplitude 
of £+ is higher for deeper surveys. Fig.^ also highlights a couple of degeneracies: for 
(z) = 1.0, now with f2 m = 0.38 is shown. Note the difficulty in distinguishing this curve 
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Fig. 1. Normalised source redshift distributions, p(z), for (z) = 0.5 (solid curve) and 
(z) = 1.0 (dashed curve) using the parameterisation of Smail et al. (1995). 

from that of the fiducial cosmology with a higher mean redshift, (z) = 1.2. The well- 
known degeneracy between f2 m and ag is also demonstrated: we show that for (z) = 1.0, 
taking for instance Q m = 0.4 and ag = 0.78 results in a curve which would be difficult 
to distinguish from that of the fiducial cosmology in practice. 



3. Estimators 

Let us now consider how estimates of the correlation functions can be obtained in practice 
using the distorted images of galaxies. We begin by considering the case where the only 
mechanism giving rise to correlations in galaxy ellipticities is gravitational lensing. 

In the limit of weak lensing, the complex ellipticity of a galaxy at position Bi is 
related to its intrinsic ellipticity e% and to the gravitational shear 7(0j) through 

e i = el+ 7 (0 i ). (14) 

We use the definition of ellipticity from Bonnet & Mellier (1995) where for elliptical 
isophotes with axis ratio b/a < 1, |e| = (I — b / a) (1 + b / a)^ 1 . Further, under this definition 
(e) ~ 7 if the unlensed source population is randomly oriented. For simplicity we will use 
the notation ji = 7(#i). Using ( f[4"| ) it follows that the expectation value 

(e l£ *) = ((4 + 7i ) (e s /+ 7 *)) = (7,7;) + (& s j) , (15) 
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Fig. 2. Dependence of £+ on angular scale 9 for (z) = 0.8 (dotted curve), (z) = 1.0 
(solid curve) and (z) = 1.2 (short dashed curve). The cosmology for these curves is as 
described in the text, with £l m = 0.3. The solid boxes correspond to (z) = 1.0, now with 
£! m = 0.38. The crosses represent (z) = 1.0, fi m = 0.4 and cr 8 = 0.78. 



and in the absence of intrinsic correlations the final term vanishes. In this case, any 
correlation between observed galaxy ellipticities results from weak lensing. 

As discussed in Schneider et al. (2002b), the correlation function can be estimated in 
bins of width A$, centred on and a weight function can be defined such that A$(<p) = 1 
for ■& — Ai)/2 < (p < {) + Ai3/2 and is zero otherwise. An estimator for the correlation 
function £+($) is given by 

um _ e, ?m a,. !«, -. llM1 NpW ^ WiWMei _ ei ^ (16) 

P ij 

where Wj and Wj are weights (depending on the reliability of the ellipticity measurement, 
for example) and iVp($) is the effective number of pairs in that bin. In the absence of 
intrinsic correlations, since (e^e*) = — it follows that is an unbiased 

estimator of the shear correlation function of the lensing signal £ + ($) . 

However, the existence of intrinsic alignments would imply that (e s ie s j) ^ 0. The esti- 
mator £ + (■$) is no longer an unbiased estimator of £ + (t?) : rather it includes a contribution 
from correlated intrinsic (source) ellipticities, for which the two-point correlation func- 
tion will be denoted by £+*($). In order to minimise the impact of intrinsic correlations 
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on the interpretation of the weak lensing signal, we need an estimate of the amplitude of 
j and a scheme to give less weight to galaxy pairs where the intrinsic correlation 
is likely to be high. 

4. Intrinsic alignments 

Correlations between the intrinsic ellipticities of neighbouring galaxies are expected to 
arise during their formation, due to similarities in the tidal gravitational field in which 
they form. As noted in the introduction, a number of numerical and analytical stud- 
ies have sought to quantify the amplitude of this correlation as a function of angular 
separation and redshift. For our purposes, we need a toy model to roughly predict the 
amplitude of intrinsic correlations expected. Since there is disagreement on the dominant 
mechanism responsible for intrinsic alignments, and on its subsequent enhancement or 
"washing out" , use of a detailed model is not warranted here. 

Let us consider a close pair of galaxies with comoving distances Wi w Wj w w and 
angular separation 9. Their comoving separation r is given by 

r 2 = (A W ) 2 + / 2 H(fl) 2 , (17) 

where Aw w Wi — Wj . 

Denoting the (3-D) two-point density correlation function by £a(r), one might expect 
that the two-point correlations in intrinsic cllipticity cx £ 2 (r). Indeed, this dependence 
corresponds to the asymptotic behaviour at large separations derived by CrOl, whose 
normalisation we later adopt. The two-point density correlation function is the Fourier 
transform of the power spectrum, so the power spectrum used to predict the lensing 
correlation functions described in Sect. 2 could be used. However, since CrOl focus on a 
£g(r) cx 1/r behaviour, and because of its analytic simplicity, we also use this form. 

Next, this 3-D intrinsic correlation function has to be projected into an angular intrin- 
sic ellipticity correlation function, integrating over the distance distribution of galaxies 
and taking into account the (observed) galaxy two-point correlation function £ gg (r) . Note 
that we do not consider the individual tangential and cross components of the correlation 
functions of the intrinsic ellipticities. Rather, we consider the sum of these quantities, 
which corresponds to £ + in the notation of Schneider et al. (2002a; note that the notation 
of CrOl differs here). We obtain 

«nt/m J dw lP ( Wl ) J dw 2 p(w 2 )[l +g gg (r)] [&(r)] 2 

fd Wl p( Wl )fdw 2 p(w 2 )[l + Z eg (r)} ' 1 
where p(w)dw is the galaxy distance distribution, which we take to be of the form given 
in Sect. 2. 

Given that the scale over which intrinsic alignments operate is much smaller than the 
distances to the galaxies, the integral above can be recast as 
,int^ JdwpHw) jd(Aw) [1+ W)] Mr)] 2 

4+ U 1 + fdwp2(w) Jd(Aw) [£ gg (r)] ' 1 > 
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The two-point galaxy correlation function is taken to be £ gg (r) = (r/4.3 /i _1 Mpc) -1 ' 8 , 
constant in comoving separation. Evolution is ignored for both the galaxy and density 
correlation functions, in view of the weak evolution that is seen in galaxy samples out to 
z <~ 3 (e.g. Porciani & Giavalisco 2002). The normalisation of is set at an angular 

scale of 10' and for (z) = 1, using the value obtained by CrOl; Mackey et al. (2002) 
point out that on very small spatial scales (<~ l/i _1 Mpc) non-linear clustering effects 
may erase any initial alignments, so we perform the normalisation well above this scale. 
We end up with a simple model for * (9) , which depends upon the redshift distribution 
of sources. 

5. Inclusion of a redshift dependent weighting factor 

An important difference between ellipticity correlations arising from lcnsing and those due 
to intrinsic alignment is that the latter only has a significant amplitude over a relatively 
small range in redshift. For a given galaxy pair, the lensing correlations depend on the 
integrated effect of the matter along the line of sight out to the redshift of the closer 
galaxy. 

In order to minimise the intrinsic signal, giving less weight to galaxy pairs where 
the redshift difference is small seems like a plausible step. How effective would such a 
weighting be, how would the lensing correlation be changed and how would the noise on 
our estimate of the correlation function increase? 

Comparison of the observed correlation functions with the foregoing expressions re- 
quires one to adopt a source redshift distribution; in practice, what will become routinely 
available are photometric redshift estimates. The true redshift of galaxy i will be denoted 
by Zi and the photometric redshift estimate by z i} with equivalent notation for distances 
(i.e. Wi and Wi). The probability density to have sources at comoving distance Wi and to 
have photometric distance estimates Wi will be denoted by p(wi) and p(wi) respectively, 
and p(wi,Wj) is the joint probability to have sources at Wi and at Wj . For conditional 
probabilities, p(u>i\wi) represents the probability that the true source distance is Wi given 
that the photometric distance estimate is Similarly, p{wi\wi) is the probability that 
a photometric distance estimate Wi will be obtained, given that the true distance is Wi. 

5.1. Estimator 

A practical estimator for the lensing correlation function should minimise the relative 
signal from intrinsic alignments: 




JVp(O) = ^WiWjZij&t (\6i - 6,1) ; (20) 



Zij = Z(zi, Zj) is a weighting factor which could be of the form 



Z(z{, zj) = 1 — exp 



r (Az) 2 i 



(21) 
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where Az is the difference between the photometric redshift estimates of the galaxies in 
a pair, and az is a factor that controls the extent of the region in Az space over which 
galaxy pairs are down-weighted. This would be chosen to reflect photometric uncertainty, 
which is typically much larger than the range over which intrinsic correlations operate. 
We will neglect the weights W in what follows (since these would be determined from 
the data in practice), and see how this estimator performs in Sect. 6. 

5.2. Leasing correlation function 

We first consider what effect a redshift dependent weighting factor has on the lensing 
induced component of the correlation function, £ + (9) . The addition of such a weighting 
factor modifies the effective source distance distribution, and a distance correlation is 
introduced. To make this clear, consider the shear correlation function for particular 
source distances Wi and wj i.e. £+(8;wi,Wj), which would correspond to sources located 
on two sheets (Schneider et al 2002a); 

9Hfo 2 dw 



Z + (9;Wi,Wj) , 

4c 4 J a z (w) 

f dll (I \ 

xR(w,Wi)R(w, Wj ) J j—P s [j^> w J J o(^): (22) 

where the upper limit Wij in the outer integral is the minimum of Wi and Wj. What 
is important to note is the dependence on the product R(w,Wi) R(w,Wj). When ( [22| ) is 
averaged over the source distance distribution, Schneider et al. (2002a) obtain 

Me) = 9J PC$k I ^«(755.-)««)W-.«.)i«(-.«.» . (23) 

where the angular brackets in (R(w, Wi) R(w, Wj)) denote averaging over the probability 
distribution p(wi,Wj). If p(wi,Wj) is the probability density to have sources at comov- 
ing distances iu, and Wj, in the absence of intrinsic source correlations one has that 
p(uii,Wj) = p(wi)p(wj). Then, ((R(w,Wi)R(w,Wj)) = W 2 (w), as in in <^). The addition 
of the weighting factor Z(wi, wj) changes the probability density of the redshifts of pairs 
included in the estimator (^p| ) thus 

, \ So" d ^ x Jo" y (OA\ 

P\W z ,Wj) — WH _ „ WH _ WH WH , ^Z4J 

Jo dw * Jo dw i z ( w ^ w j)Jo dw i x j dw jy 

where X = p(wi) p(wi\wi) and y = p(u)j) p(Wj\wj), the inclusion of the denomina- 
tor ensuring that the probability distribution is normalised. Using p(wi\wi) p{wi) = 
p(wi\wi) p(u>i) (and similarly for galaxy j) we have that the expectation value of 

, . v . f w " dwi R(w,Wi) f WH dwj R(w,Wj) L H dm,* X' f™ H dw, ; y' , 

(Riw^R^Wj)) = ^ Xh a % A 3 f U' - t r WH \_ Z — ,(25) 

Jo dw * Jo dw i Jo dwi X Jo dw i y 

where X' = p(wi) p(wi\wi) and y' = p(wj) p(wj\wj)Z(wi,Wj). Substitution of ( [25] ) 

into (|23| ) gives the lensing two-point correlation (9) . This depends upon the accuracy 

of photometric redshift estimates, and upon "filtering" using the weighting factor Z. 
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5.3. Intrinsic correlation function 

Now consider how accounting for photometric redshifts and applying the same weighting 
factor Z(wi,Wj) to the prescription for the intrinsic correlation function adjusts the 
intrinsic correlation function: 



where X' and y are defined as below fl2q). The £gg( r ) term in the denominator can 
safely be ignored, leaving us with the same denominator as (^5j). We note that the 
galaxy correlation function £ gg should in principle also enter the expressions for the 
shear; however, as shown in Schneider et al. (2002a), this modifies the resulting cosmic 
shear signal only by a very small fraction and has therefore been ignored throughout. 



The error in each of our photometric redshift estimates is assumed to be a Gaussian of 
dispersion cr p hot = 0.1 in z, typical of the accuracy obtained using codes such as hyperz, 
which adopts a standard spectral energy distribution fitting procedure (Bolzonclla ct 
al. 2000), provided that photometric information is available in a sufficient number of 
wavelength bands. The ratio of the intrinsic signal to the total (lensing plus intrinsic) 
signal for the weighting factor given by ( pl| ) and for various widths of the weighting 
function, az, is shown in Fig.[| Two redshift distributions are shown, one with (z) = 1.0 
and the other with (z) = 0.5. Note that the contamination from intrinsic alignments is 
much more dramatic in the lower redshift case, but that this can be largely suppressed 
by down- weighting pairs of galaxies with similar photometric redshifts. The higher mean 
redshift case is perhaps more typical of surveys suitable for cosmic shear studies. Although 
the fractional contribution from the intrinsic alignment signal is low, it can be reduced. 

By down-weighting pairs depending on their photometric redshift estimates, the num- 
ber of pairs is also effectively reduced and hence the noise is increased. For the higher 
(lower) redshift case, the solid (dashed) curve of Fig.|| shows the fraction of pairs re- 
maining (J 7 ) as a function of the filter width crz- Although the noise is increased, the 
expectation value (R(w,Wi) R(w,Wj)) is rather insensitive to <Jz and hence the change 
in £+(#) is also small (see Fig.||). 

7. Discussion and Conclusions 

Recent studies of intrinsic ellipticity correlations have shown that the intrinsic signal 
can be comparable to (or exceed) the lensing signal in the case of shallow cosmic shear 
surveys (Jing 2002). We investigated to what extent the contribution of galaxy pairs 
which are likely to have larger intrinsic correlations can be suppressed, using photometric 
redshift information and a simple redshift dependent weighting factor. Such a process also 




(26) 



6. Results 
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Fig. 3. The ratio of the intrinsic correlation function £™* to the combined lensing and 
intrinsic correlation function (£ + + as a function of angular scale 8. For all cases, 
the photometric redshift error is assumed to be Gaussian distributed, with dispersion 
Cphot = 0.1. The lower set of three curves are for (z) = 1.0, with the curve marked "A" 
corresponding to the case where there is no redshift dependent weighting factor, Z. The 
curve below this is for oz = 0.1 and the lowermost curve is for oz = 0.4. The upper 
set of three curves are for a lower mean redshift (z) = 0.5, with the curve marked "B" 
being the case where there is no weighting factor. As for the higher mean redshift case, 
the other two curves correspond to oz =0.1 and oz = 0.4 respectively. 



results in a reduction in the effective number of pairs and hence an increase in the noise, 
which scales roughly as 1 / y/N p . The width of the weighting factor, parameterised by 
<7z, controls the fraction of pairs remaining. The photometric accuracy determines the 
precision with which galaxies at similar redshifts can be identified, hence the residual 
contribution from intrinsic alignments. 

The reduction of the effective number of pairs is not the only relevant consideration 
for the noise in the determination of the two-point correlation function. Depending on the 
survey geometry, cosmic variance can become the dominant contribution to the covariance 
of the correlation function for larger angular separations. In particular, for a compact 
survey geometry, cosmic variance of the shear field dominates over shot noise for angular 
scales larger than a few arcminutes (e.g., Schneider et al. 2002b). Hence, for those scales 
the reduction of the effective number of pairs becomes of little relevance. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Fig. 4. The effective fraction of pairs remaining (J 7 ) as a function of oz (filter width) 
for (z) = 1.0 (solid curve) and (z) — 0.5 (dashed curve). 

Intrinsic alignments have also been invoked as a possible source of the so-called 
B-mode contributions seen in cosmic shear surveys (van Waerbeke et al. 2001, 2002; 
Hoekstra et al. 2002b). Such B-modes are not expected from lensing effects and thus 
are usually interpreted as a remaining systcmatics; whereas galaxy correlations can in 
principle generate a B-mode contribution (Schneider et al. 2002a), its amplitude is very 
small. If the B-mode is due to an intrinsic alignment of galaxies, then the use of the 
redshift filter as discussed here should strongly suppress its relative contribution. 

Choosing to apply such a weighting factor will be governed by the details of the survey. 
Such considerations include that the impact of intrinsic alignments is more dramatic for 
surveys of lower mean redshift or depth. Another important factor is the accuracy of 
photometric redshift estimates, which depends upon the combination and characteristics 
of filters used for the survey - see e.g. Wolf et al. (2001) for a discussion of photometric 
redshift estimate performance in the context of filter sets. Further, the size of the error 
bars associated with the practical determination of the two-point correlation function will 
also play a part in the decision. For deep surveys, with our assumed photometric errors, 
in order to decrease the fractional contribution of ^Y 1 * by a few percent, a reduction of 
several tens of percent in the number of effective pairs may result. However, removal of the 
intrinsic alignment systematic becomes increasingly important as the experimental error 
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Fig. 5. The lensing correlation function £+(9) as a function of 9, for (z) = 1 (upper 
curves) and for (z) = 0.5 (lower curves). Within each set of three curves, the uppermost 
curve is for no weighting, the middle curve is for oz =0.1 and the lowermost curve is 
for o z = 0.4. 

bars in surveys become smaller. In shallower surveys, the contamination from intrinsic 
alignments is much more pronounced, and in this case down- weighting close galaxy pairs 
is more effective. Given real data, one might consider a more elaborate weighting scheme 
that is a function of not only the difference in photometric redshifts, but also giving more 
weight to higher rcdshift pairs. Furthermore, the weight factor Zij need not be chosen 
as a function of the estimated photometric redshift only, but can be constructed such 
that it depends on the full rcdshift probability distribution as estimated by photometric 
redshift techniques. For example, a natural choice would be to have Zij depending on 
the probability that the two galaxies i,j lie within a narrow redshift interval Az, as 
estimated from their individual z-distributions. The calculation of the expectation value 
of the corresponding shear estimator is then slightly more complicated, but can be done 
for each survey at hand. 

Throughout we have assumed that photometric redshift estimates will be available 
for each of the source galaxies, and that the dispersion in these estimates is independent 
of magnitude or redshift. In practice, photometric rcdshift estimates are reliable down 
to a certain limiting magnitude, depending on the survey characteristics and the relative 
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depths of the filters used. Particularly in a deep survey, this magnitude limit may exclude 
many of the fainter and smaller galaxies that are likely to be in the high-redshift tail of 
the redshift distribution, and to make a substantial contribution to the lensing signal. 
Requiring photometric redshift estimates could in fact lead to an increase in the Poisson 
noise of such a survey by reducing the usable number of pairs beyond the effect shown 
in Fig. 4, and in addition the reduction in mean redshift could increase the relative 
contribution from intrinsic alignments. One possibility to counteract such affects might 
be to include fainter galaxies where no photometric redshift estimates are available, 
assuming that these are at higher redshifts where intrinsic alignment is less of an issue. 

The use of cosmic shear surveys as a tool to place constraints on the 3-D power spec- 
trum, so-called power spectrum tomography, would rely on the availability of photometric 
redshift estimates for the galaxies involved in the analysis (e.g. Hu 1999; 2002). Croft 
& Metzler (2000) demonstrated that the intrinsic alignment signal is more severe when 
correlations within a narrow slice in source galaxy redshift space are considered. The ap- 
plication of a redshift dependent filter, similar to the one we have described here, would 
greatly suppress the intrinsic alignment systematic. Instead of slicing galaxies in redshift, 
and calculating the correlation function for galaxies in the same redshift bin - which both 
increases the relative importance of intrinsic alignments and reduces the total number of 
pairs - one should correlate galaxies from different redshift bins i,j, where the bin width 
would be chosen to be of the same order as the redshift uncertainty. This would yield a 
set £ij (9) of correlation functions for which the intrinsic signal is strongly suppressed for 
i 7^ j. When comparing these correlation functions with predictions from cosmological 
models, their covariance must be taken into account, which may turn out to be fairly 
complicated (see Schneider et al. 2002b for the case with no redshift information). This 
redshift slicing is most straightforwardly done with the correlation function, although for 
other estimators of the two-point cosmic shear statistics, such as the shear dispersion or 
the aperture mass, similar pair redshift-dependent estimators are easily constructed. 

Most of what has been said above also applies to higher-order cosmic shear statistics. 
It is well known that the three-point statistics, i.e. the skewness, contains very useful 
cosmological information (e.g. van Waerbeke et al. 1999). As true for the 2-point statistics, 
the three-point correlation function is the quantity which is most easily derived from a 
cosmic shear survey. Bernardeau et al. (2002a) have constructed a statistics based on 
the measured three-point correlation function (see Schneider & Lombardi 2002 for the 
classification of three-point shear correlators), and Bernardeau et al. (2002b) obtained a 
significant detection in the VIRMOS-DESCART survey data. It is unclear whether, and 
by how much, intrinsic cllipticity correlations affect these measurements. Furthermore, 
up to now no B-mode estimator in the three-point function has been devised which may 
indicate the presence of intrinsic alignment effects. Redshift slicing can of course also be 
done for the three-point function which may be the only way to measure these lensing 
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statistics without the potential influence of intrinsic effects; one needs to suppress those 
triplets of galaxies where the probability of all three being at the same redshift is not 
negligibly small (i.e., the estimator is unaffected by intrinsic effects if two of the three 
are at the same redshift). 

Of course, the intrinsic alignment signal that we seek to suppress when performing 
a lensing analysis is also interesting in its own right, since it places constraints on the 
formation and evolution of galaxies. Suitable catalogues of galaxies to quantify this signal, 
and its evolution with redshift, will be a natural by-product of large cosmic shear surveys 
with photometric redshift information. For example, using a filter for the determination 
of the ellipticity correlation as in (|2(]) with Z[ a = 1 — -Hjj will make this estimate 
dominated by the intrinsic alignments, and could thus be used as a first estimate of it. 
More sophisticated techniques would include the consideration of the resulting correlation 
function in dependence of the width of the weighting function, and extracting the lensing 
and intrinsic signal from this functional form. 

If the intrinsic alignment of galaxies indeed occurs, and in particular if the intrinsic 
effect is as strong as suggested by some models (e.g. Jing 2002), a deep multi-colour 
wide-field survey with a broad range of filters will be necessary, and rewarding, to remove 
this systematic from cosmic shear measurements. Most likely, the near-IR imaging will 
present the bottleneck, limiting the magnitude - and thus the effective number density - 
of galaxies that can be used for cosmic shear. A wide-field near-IR camera in space, such 
as the PRIME satellite mission, would be an ideal supplement to the planned extensive 
ground-based optical cosmic shear surveys. 

After we had completed this paper, we became aware of work by Heymans & Heavens 
(2002), also addressing how redshift information can be used to reduce the contamination 
from intrinsic alignments. They apply their technique to estimate the reduction of the 
intrinsic signal in several surveys, including the SDSS photometric and spectroscopic 
samples. 
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